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Abstract. In order to assess the potential of future mi- 
crowave anisotropy space experiments for detecting clus- 
ters by their Sunyaev-Zeldovich (SZ) thermal effect, we 
have simulated maps of the large scale distribution of their 
Compton parameter y and of the temperature anisotropy 
AT/T induced by their proper motion. Our model is based 
on a predicted distribution of clusters per unit redshift and 
flux density using a Press-Schecter approach (De Luca et 
al. 1995). 

These maps were used to create simulated microwave 
sky by adding them to the microwave contributions of the 
emissions of our Galaxy (free-free, dust and synchrotron) 
and the primary Cosmic Microwave Background (CMB) 
anisotropies (corresponding to a COBE-normalized stan- 
dard Cold Dark Model scenario). In order to simulate 
measurements representative of what current technology 
should achieve, "observations" were performed accord- 
ing to the instrumental characteristics (number of spec- 
tral bands, angular resolutions and detector sensitivity) of 
the COBRAS/SAMBA space mission. These observations 
were separated into physical components by an extension 
of the Wiener filtering theory (Bouchet et al. 1996). We 
then analyzed the resulting y and AT/T maps which now 
include both the primary anisotropies and those superim- 
posed due to cluster motions. A cluster list was obtained 
from the recovered y maps, and their profiles compared 
with the input ones. Even for low y-values, the input and 
output profiles show good agreement, most notably in the 
outer parts of the profile where values as low as y ~ 3. 10 -7 
are properly mapped. We also construct and optimize a 
spatial filter which is used to derive the accuracy on the 
measurement of the radial peculiar velocity of a detected 
cluster. We derive the accuracy of the mapping of the very 
large scale cosmic velocity field obtained from such mea- 
surements. 
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1. Introduction 

Future space experiments should allow to built a large 
statistically homogeneous sample of clusters of galaxies 
detected by the spectral signature of CMB photons scater- 
ring off the free electrons of the hot intra-cluster gas. This 
would be of great cosmological interest, allowing to con- 
strain the cosmological scenarios of large scale structure 
formation and evolution, as well as the gas history (Co- 
lafrancesco & Vittoriol994, Barbosa et al. 1996). 

The measurement of the peculiar velocity of clusters 
of galaxies could be an important tool to study the large 
scale velocity field of the universe. This in turn provides a 
unique opportunity for probing the underlying mass distri- 
bution. Thus one can probe fairly directly the primordial 
spectrum and further constrain the various cosmological 
models. The direct determination of the peculiar veloc- 
ity (by independent redshift and distance determination) 
is observationally difficult and time-consuming. Still, sev- 
eral groups have succeeded in measuring the bulk, volume- 
averaged, peculiar velocity field, in our neighborhood, on 
scales r = 10 h^ 1 to 60 hr 1 Mpc (Aaronson et al. 1986, 
Collins et al. 1986, Dressier et al. 1987, Strauss & Willick 
1995 (for a recent review)). 

Another method is nevertheless possible. It relies on 
combining the informations from the Sunyaev-Zeldovich 
(SZ) effects (thermal and kinetic) which create secondary 
temperature fluctuations (Zcldovich & Sunyaev 1969, Sun- 
yaev & Zeldovich 1972, 1980). The combination of the 
measurements of the two effects give a fairly direct de- 
termination of the peculiar radial velocity of clusters of 
galaxies (Sect. 2.2, Eq. 3), although a good estimate of 
the peculiar velocity requires a precise measurement of 
both SZ effects and of the gas temperature. A crucial ad- 
vantage of such a combination is that the SZ effects do not 
decrease in brightness with distance. Therefore measure- 
ments of both effects for a large number of clusters could 
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give precise estimates of the large scale velocity field dis- 
tribution. 

The temperature fluctuations generated by the kinetic 
SZ effect have the same spectral signature than the pri- 
mordial anisotropies of the CMB. For cluster velocity de- 
terminations, the primordial temperature fluctuations of 
the CMB act as a contaminating source of the SZ kinetic 
effect. In addition, one must also take into account all 
other sources of contamination that might spoil the mea- 
surement, from resolved sources (other clusters, galaxies 
...), unresolved ones (galactic synchrotron, free- free emis- 
sion ...) or instrumental noise. 

Some measurements of the SZ thermal effect, on known 
clusters, have already been made in the wavelength range 
where the SZ effect creates a decrement in the CMB in- 
tensity (A > 1.3 mm) (see Rephaeli 1995 and references 
therein). 

In the context of the feasibility study of future space 
projects such as the COBRAS/SAMBA missiontl, of the 
European Space Agency, dedicated to the CMB observa- 
tions together with other cosmological targets (SZ detec- 
tion in clusters, primordial galaxies, ...), a complete simu- 
lation of the astrophysical processes, the expected instru- 
mental characteristics of the satellite and of the separation 
of these processes was performed (Bouchet et al. 1996). 
This simulation was used as a tool to constrain and quan- 
tify the capabilities of the satellite. In this paper, we focus 
on the capabilities of such a mission for the detection of 
the SZ thermal effect of clusters of galaxies, the imaging of 
the clusters and finally the measurement of their peculiar 
velocities. 

Section 2 introduces briefly the formalism of the ther- 
mal and kinetic SZ effects, while we present in section 
3 the method we used to generate SZ effect maps. Sec- 
tion 4 assesses the capabilities of a space experiment like 
COBRAS/SAMBA to detect clusters, and we evaluate, in 
section 5, the accuracy of the peculiar velocity determina- 
tion of individual clusters. We finally derive the expected 
accuracy on the rms dispersion of the large scale velocity 
field in section 6. Section 7 summarizes the main results 
and conclusions of the paper. 



eter y, which depends only on the cluster's electronic tem- 
perature and density (T e , n e ): 



m e c* 



T e (l)n e (l)dl 



where k is the Boltzmann constant, o~t the Thomson cross 
section, m e c 2 the electron rest mass energy and I is the 
distance along the line of sight. When the intra-cluster 
gas is isothermal (T e (l) — T e ~ const), y is expressed as a 
function of the optical depth t (t = ctt f n e (l) dl): 
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The relative monochromatic intensity variation of the 
CMB due to the SZ thermal effect is given by 



= V ■ f{x) 



where x is the adimcntionnal frequency x = hv /kTcMB 
(h denotes the Planck constant, Tqmb the CMB tempera- 
ture, and v the frequency) , /„ is the intensity of the CMB 
(black body emission) and f{x) is the spectral form factor 
given by: 



(e x - 1) 
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2.2. Kinetic Sunyaev-Zeldovich effect 

If a cluster has a radial peculiar velocity v r , another rel- 
ative intensity variation of the CMB due to the Doppler 
first-order effect is added. It is given by: 



^ = ~T-a{x) = I -jr ) -ai.r). 



(AT 



(2) 



SZ 



where a(x) is the spectral form factor for the kinetic effect, 
given by: 

a{x) = x- 



1 



2. Sunyaev-Zeldovich thermal and kinetic effects 

2.1. Thermal Sunyaev-Zeldovich effect 

The thermal SZ effect is the inverse Compton scattering 
of CMB photons by free electrons in the hot intra-cluster 
medium. Since the number of photons is conserved, their 
spectrum is just shifted on average to higher frequencies. 
This effect is characterized by the comptonization param- 



and r is the optical depth. The effect is positive for clusters 
moving towards the observer (i.e. with negative velocities). 
The intensity fluctuation induced by the SZ kinetic effect 
has the same spectral shape as the primordial ones (equiv- 
alent to a temperature fluctuation). 

The temperature variation due to the kinetic SZ effect 
can be written as a function of the y parameter (defined 
in Eq. |l|). The combination of the two SZ effects gives the 
expression of the radial velocity of the cluster: 



1 For a complete description of the COBRAS/SAMBA 
project, see ESA document D/SCI(96)3. 
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At high frequencies, the expressions given above are not 
accurate enough and relativistic calculations of the ther- 
mal effect are needed in most cases (Rephaeli 1995 and 
references therein) . A relativistic treatment introduces dif- 
ferences in the relative intensity variation and a shift of 
the crossover frequency. The corrections depend on both 
the temperature of the intracluster medium and the fre- 
quency. In this paper, we have restricted our study to the 
nonrclativistic treatment as a "text book" case in order to 
adress the capabilities of a space mission in measuring the 
SZ effect on clusters of galaxies. The corrections for the 
relativistic case will have to be taken into account when 
dealing with real data. 

3. Simulations 

3.1. Cluster model 

We model the y profile of a single resolved cluster ac- 
cording to a King model (King 1966). The main physical 
characteristics of the cluster are n e (R) and T e (R), respec- 
tively the electronic density and temperature distributions 
given as functions of the distance R to the center of the 
cluster. Hereafter, we use a hydrostatic isothermal model 
with a spherical geometry (Cavaliere & Fusco Femiano 
1978, Birkinshaw, Hugues & Arnaud 1991). For the den- 
sity distribution, this gives: 



n e (R) = n e0 



R c 



3/3 



where n e o is the central electronic density, R c is the core 
radius and j3 is a parameter of the model which represents 
the ratio of the kinetic energy per unit mass in the galaxies 
to the one in the gas. We take (3 = 2/3, as indicated 
by both numerical simulations (Evrard 1990) and X-ray 
surface brightness profiles (Jones & Forman 1984, Edge & 
Stewart 1991). The hypothesis of isothermality is rather 
well confirmed by the X-rays observations of the ASCA 
satellite (Mushotzky 1994). 

In this specific model for the gas distribution, the in- 
tegrated y profile of the cluster is given by: 



V oc 



1+ ^ 



1 -1/2 



9 and 9 C being respectively the angular distance to the 
center of the cluster and core radius. The Full Width at 
Half Maximum (FWHM) of the profile and the core radius 
of the cluster are related by FWHM = 3.5 x 9 C . 

For the time evolution of the temperature T e and core 
radius R c , we use the model of Bartlett and Silk (1994) in 
which the main parameters of the cluster evolve as sug- 
gested by the self-similarity arguments (Kaiser 1986), with 
a parameter standing for the negative evolution of the 
number of clusters indicated by the X-ray observations. 



The normalization of the temperature and core radius is 
made using the A665 ASCA data (Yamashita 1994). 

3.2. The maps 

We simulate 500 x 500 pixels (12.5° x 12.5°) maps of the 
clusters, in terms of their y parameters for the thermal 
SZ effect and AT/T maps of the same size for the kinetic 
SZ effect. The distribution of clusters per unit of redshift, 
solid angle and flux density interval was computed (De 
Luca ct al. 1995) using a Press-Schechter mass function 
(Press & Schechter 1974) normalized to X-ray and optical 
data (Bahcall & Cen 1993). The total number of sources 
in each map is drawn using a Poisson distribution, the 
position of each cluster being also assigned at random, i.e. 
the effect of spatial correlations are not taken into account. 
The counts in De Luca et al. (1995) are obtained using a 
formalism similar to the one adopted in Bartlett & Silk 
(1994) taking into account the negative evolution in time 
suggested by X-ray observations. 

A discrimination between resolved and unresolved 
clusters is made according to their spatial extent. We 
give the gas distribution within clusters a finite extent 
with a maximum radius R ma x- In fact, the hydrodynamic 
isothermal model is in good agreement with X-ray obser- 
vations over 5 core radii (Markevitch et al. 1992), other 
observations (Henriksen & Mushotzky 1985) show that 
in some cases the gas in clusters is still seen up to 10 
core radii. This seems to be also the case in the ASCA 
data (Mushotzky 1994). These sizes are the ones estimated 
from the X-ray observations. Since the X-ray emisssion is 
proportional to n^, the observations are strongly biaised 
towards the center of the clusters and tend to minimize the 
extent of the clusters. Recent simulations of cluster forma- 
tion (Evrard et al. 1996) show that the virialization limit 
for the clusters takes place at 2 to 3 Mpc from the center. 
For typical values of the core radius (0.1 to 0.2 Mpc) the 
dynamical extent of the clusters is thus of the order of 10 
to a few tens of core radii. Because the SZ effect is propor- 
tional to the n e , it is therefore very sensitive to the outer 
parts of the density profile where most of the mass rests. 
In order to investigate the case where the clusters extend 
beyond the estimations from X-rays, we assume, for the 
simulated maps, that the maximum extent of a cluster is 
2W C with a profile going down to zero for 9 > 1Q9 C . In the 
maps, the unresolved sources are thus the ones for which 
2i)9 c < 1.5', where 1.5' is the simulation pixel size. In that 
case, we assign to the corresponding pixel an integrated 
y parameter, noted Y (see Eq. 4). When the sources are 
spatially resolved (2Q9 C > 1.5'), we compute their y pro- 
file using the above /3-model for the gas (Sect 3.1 and 
appendix A). 

To generate the AT/T maps arising from the SZ ki- 
netic effect, the radial peculiar velocity v r of each cluster is 
drawn randomly from an assumed Gaussian velocity dis- 
tribution with standard deviation today er = 400 km/s 
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(Faber et al. 1993). The time evolution of the standard 
deviation is followed according to linear perturbation the- 
ory. We also assume no correlations in the velocity distri- 
bution. Using the same procedure than above concerning 
the source extents, we compute AT/T profiles for the re- 
solved clusters (Appendix A). 

We found that this simple model based on De Luca et al. 
(1995) modeling for the SZ source counts turns out to be in 
rather good agreement with Bond & Myers' (1996) more 
sophisticated simulations, based on the "peak patch" al- 
gorithm, at least for statistical quantities such as the rms 
values of both the y parameter and AT/T. This justifies 
the approximations made concerning correlations. 

4. Detection and mapping of distant clusters 

It is now possible, due to recent technological improve- 
ments, to design a new generation of satellites dedicated 
to CMB observations at small angular scales (few arcmin- 
utes to few degrees) with very high sensitivities, allow- 
ing the detection of temperature fluctuations at a level of 
AT/T ~ 10~ 6 . Measuring primordial temperature fluc- 
tuations at this level of accuracy requires the ability to 
detect secondary fluctuations such as the ones associated 
with the SZ effect. Therefore one of the byproducts of such 
a survey of the CMB will certainly be a new catalogue of 
clusters of galaxies. 

The physics of cluster formation and the history of gas 
virialization are still almost unknown. The SZ effect mea- 
surement provides a new method of observing clusters 
which is potentially more powerful than X-ray observa- 
tions for the search of clusters at high redshifts. 
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Fig. 1. Ratio of the X brightness to the intensity of the SZ 
thermal effect as a function of the redshift, for three tempera- 
tures of the cluster. 

Figure 1 gives the ratio of the X-ray brightness of a galaxy 
cluster to the intensity of its SZ thermal effect, as a func- 



tion of the redshift. It shows that at z = 1 the X-ray 
brightness drops by a factor greater than 30 compared to 
the SZ intensity, indicating that SZ observations should 
be very powerful to detect high z clusters. 

In the context of CMB observations at small angu- 
lar scales with high sensitivities, the problem of the de- 
tectability of clusters among the other astrophysical com- 
ponents is essential to evaluate how many clusters could 
be found in a sky survey, especially at large distances. 
Fortunately, the SZ thermal effect, which has a specific 
spectral signature (positive and negative intensity peaks 
at respectively 0.85 and 2 mm with a zero value at 1.38 
mm), is rather easy to identify in the framework of a multi- 
frequency experiment covering the range 30 to 800 GHz, 
with high sensitivities and good angular resolution such 
as COBRAS/SAMBA. Using the separation of processes 
(primordial CMB, clusters, free-free, ...) which is based on 
the projection of the "observed" signal in different wave- 
bands (Bouchet et al. 1996), the sky simulation provides 
us with a recovered y map for the SZ thermal effect to- 
gether with a AT/T map. These recovered maps are used 
to characterize the number of clusters that could be de- 
tected and furthermore they show the ability of the survey 
to image the clusters. 

In order to find the detected clusters in the recovered y 
maps, we use the following algorithm. Above a threshold 
of y = 2. 10 -6 , each maximum is associated with the cen- 
tral peak of the y profile of a resolved cluster. Assuming 
spherical symmetry, we compute the radial profile of the 
comptonization parameter y, by averaging within rings of 
equal width. We also compute the integrated y parame- 
ter, Y, over the profile for each recovered cluster. Then, we 
compare these reconstructed profiles with the input ones. 
We thus produce a catalogue of detected clusters giving 
the position of the y maxima in the map, together with 
the "measured" central value yo, and the integrated pa- 
rameter Y. The method, which consists of integrating the 
signal over rings getting larger and larger, tends to over- 
estimate the signal if one does not take out a base line 
which is the average y value given by undetectable weak 
clusters. For rings of constant SR/R, both noise and signal 
decrease as 1/R where R is the distance to the center of 
the cluster; thus the signal (integrated over rings) to noise 
ratio is constant as long as the density follows a /3-profile. 

Figures 2 and 3 show both input and recovered profiles 
for respectively a strong cluster (in terms of the y param- 
eter) and a weak one, very close to the sensitivity limit 
of the COBRAS/SAMBA instruments. In both figures we 
note that the central part of the cluster suffers from the 
beam dilution, this region is therefore not recovered in 
a satisfactory way by the SZ observations. On the other 
hand, the wings of the cluster are well recovered, they are 
observable down to y ss 3. 10~ 7 over about one degree as 
indicated in Fig. 2. The comparison between the input 
and recovered profiles shows a good agreement. In fact, 
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Fig. 2. Input (solid line) and recovered (dashed) profiles for 
a cluster with high y as a function of the distance to center 
in arcminutes. The integrated y parameter is recovered with 
better than 2% 
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Fig. 3. Input (solid line) and recovered (dashed) profiles for a 
cluster 100 times weaker than typical observed clusters. Here, 
the intagrated y parameter is recovered with an accuracy of 
12%. 



the sensitivity is such that clusters, and specifically the 
wings of the profile, can be observed through the SZ ther- 
mal effect as far as the gas extends; the main limit to such 
a measurement is likely to be the confusion limit, due to 
the overlap of weak clusters in the background. The abil- 
ity to detect the wings of the profiles indicates that the SZ 
effect is a powerful tool to constrain gas accretion models 
in potential wells of clusters and measure the virialization 
radius which gives the spatial extent of a cluster. Further- 
more, more information in the statistical properties of the 
distribution of the comptonization parameter y can be ex- 
tracted and used to constrain evolution models even when 
confusion sets in. 



The X-ray emission and the SZ effect are dominated by 
different parts of the cluster; X-rays are mainly associated 
with the core of the cluster because the intensity of brem- 
strahlung emission is proportional to n 2 , whereas the SZ 
effect which is proportional to the electron column den- 
sity is dominated by contributions coming from the lower 
density regions where the path length is longer. Therefore, 
the information brought by both X-rays and the SZ effect 
are quite different and rather complementary. 
A comparison between the capabilities of X-ray and SZ 
measurements for clusters can be done by using the char- 
acteristics of the EPIC camera planned for the X-ray ob- 
servatory mission XMM and the COBRAS/SAMBA char- 
acteristics. Due to their good resolution and because they 
are biaised by the rt 2 dependence, X-ray observations are 
well adapted to map the core of the clusters. The EPIC- 
XMM instruments can map the cluster up to about 79 c at 
a sensitivity of 5<r reached in 20 hours of integration. This 
is illustrated in figure 4 for an A496-like cluster (T e o = 4 
keV at z = 0.2). The EPIC instrument will map the clus- 
ters (especially the central part) and will give the tempera- 
ture distribution, electronic density... The SZ observations 
suffer from a significant beam dilution of the core, because 
of the lack of resolution compared with the X-rays mea- 
surements, but they resolve the wings of the cluster profile 
up to the virialization limit (cf. Fig. 2). Figure 4 shows 
that the A496-like cluster is mapped up to about 17 9 c 
with SZ observation. Therefore, a combination of X-ray 
measurements together with SZ observations will strongly 
constrain cluster models. 

A high sensitivity SZ sky survey like the one that can 
be carried out by a mission like COBRAS/SAMBA will 
give the best catalogue of distant clusters and will be used 
to define follow-up observations with X-ray observations 
which require a long integration time to study a distant 
cluster. 

One can characterize the clusters in terms of their mea- 
surable parameters: angular core radius C , central value 
of the comptonization parameter yo or the integrated pa- 
rameter Y, which combines the two previous quantities. 
Assuming a King profile, this parameter is directly related 
to the mass of the gas, the redshift z and the temperature 
of the clusters, by the following expression (for a flat uni- 
verse): 



/ 



y dCl 



10 



" 4 arcmin 2 



= 0.43/1' 



( M G \ f kT e \ 
\IO 14 M ) \10 keV] 



(1 + zf 



(4) 



In terms of observable parameters, Y is given by Y = 
64yo($c/arcmin) 2 . We expect that the extraction of the 
clusters will depend mainly on the value of Y, because 
the sensitivity of the SZ profile does not decrease with 
radius as long as the gas density decreases as R~ 2 .Thus, 
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Fig. 4. Comparison between X-ray brightness (solid line) and 
the y parameter (dotted line) for a A496-like cluster (T e o = 4 
keV at z = 0.2) with XMM-EPIC and COBRAS/SAMBA res- 
olutions, respectively for the X-rays and y observations. Left 
axis: X-ray brightness in EPIC counts/s/arcmin 2 ; right axis: y 
paramter, versus angular size in arcmin. The solid horizontal 
line stands for the X-ray background, the sensitivity limit of 
the XMM instrument for an integration time of 20 hours is in- 
dicated by the dashed horizontal line whereas the dotted line 
is for the COBRAS/SAMBA sensitivity limit. 



Fig. 5. Fraction of the recovered clusters from a recovered map 
after component separation, for three sensitivities as a function 
of the integrated Y comptonization parameter. The nominal 
mission of COBRAS/SAMBA is 14 months of observations. 

dependent and could increase by a factor 3 for cosmolog- 
ical models with low f2 due to the existence of distant 
clusters (z > 0.5) and their contribution to the counts as 
in the model of Barbosa et al. 1996. 



the sensitivity required to detect a cluster is controled by 
the integrated emission and not by the peak brightness. 

We use the extraction method described above, applied 
to the simulated input and recovered maps and derive the 
completeness of the obtained catalogue of resolved clusters 
as a function of their Y. 

A systematic extraction of the clusters is done for three 
noise configurations (the nominal observing time for CO- 
BRAS/SAMBA mission is about one year). Both 48 and 
12 months of observations enable the detection of at 
least 95% of the clusters with Y > 9. 10 -4 arcmin . For 
the 12 months configuration, 70% of the clusters with Y 
> 5. 10~ 4 arcmin 2 and ya > 3. 10 -6 arc recovered. For 
such clusters, the integrated y parameter is recovered, on 
average, with an accuracy of about 30% (clusters shown 
in figures 2 and 3, the accuracies are respectively 1.4% 
and 12%). Decreasing the noise level leads obviously to 
a deeper survey as indicated in the plot (Fig. 5). If we 
increase the noise level to a noise corresponding to a 6 
months mission, even some of the strongest clusters are 
not detected. In fact, we recover only 85% of clusters for 
which Y > 10 -3 arcmin 2 , and the extraction gets worse 
for clusters with Y = 5. 10 -4 arcmin 2 . 

By using the De Luca et al. (1995) source counts 
model, we expect that a mission like COBRAS/SAMBA 
could detect about 7. 10 3 resolved clusters with Y > 
9. 10~ 4 arcmin 2 and a completeness better than 95%, and 
about 10 4 clusters with Y > 5. 10 -4 arcmin 2 and com- 
pleteness better than 68%. These numbers are very model 



5. Peculiar velocity measurement 

Given the radial peculiar velocities at known positions and 
under the assumption of a potential velocity flow around 
our galaxy, one can extract the velocity field and derive the 
bulk velocity Vr which is the average of the local velocity 
field smoothed over a window function of scale Rh~ l Mpc. 
One can get the peculiar radial velocities using the rcdshift 
surveys and some relationships, giving the distances to the 
objects, such as the Faber-Jackson (1976) or Tully-Fishcr 
(1977) relations. Several authors have measured the bulk 
velocities averaged over different volumes (a good approx- 
imation at very large scales), Dressier et al. (1987) found 
V6o = 599 ± 104 km/s, Courteau et al. (1993) measured 
V i0 = 335 ± 38 km/s and V w = 360 ± 40 km/s, Willick 
et al. (1996) computed the bulk velocity using POTENT 
they found Vqq = 20 — 300 km/s. The main difficulty from 
which this method suffers is that it requires very reliable 
distance indicators. Therefore, the relative error in deter- 
mining the peculiar velocity increases with distance. 

The determination of the peculiar velocity from SZ 
measurement is promising, since its accuracy is distance 
independent. If indeed one can detect the kinetic SZ effect 
for a number of clusters as large as 10 4 (as possible with 
the COBRAS/SAMBA mission), this will provide unique 
statistical information on the velocities. 
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5.1. Calculation method and geometrical filter 

The radial peculiar velocity of the cluster is given by com- 
bining the measurements of both kinetic and thermal SZ 
effects (Sect. 2.2, Eq. 3). The components separation takes 
advantage of the spectral signatures of the different astro- 
physical contributions (free-free, synchroton, clusters of 
galaxies, ...), which are taken into account in the sky sim- 
ulation (Bouchet et al. 1995), and of the characteristics of 
the instrument to give recovered maps of the astrophysi- 
cal processes. In particular, the AT/T map, obtained af- 
ter the components separation, includes both primordial 
CMB fluctuations and (AT/T)sz because they have the 
same spectral signature. The (AT/T)sz is thus contam- 
inated by the primordial temperature fluctuations of the 
CMB, this contamination beeing responsible for an error 
in the determination of the peculiar radial velocity of a 
cluster when one uses a method based on the combination 
of measurements of both SZ effects. The aim is to have 
the smallest error on the velocity, this requires a good 
measurement of (AT /T)sz- 

The induced temperature fluctuations due to the SZ 
kinetic effect have angular sizes smaller than one degree, 
whereas the CMB spectrum exhibits a "Doppler" peak at 
about one degree. Within this framework, a good mea- 
surement of the SZ kinetic effect (and thus the peculiar 
velocity) must be a compromize between on the one hand 
maximizing the signal by integrating over a large beam, 
and on the other hand minimizing the spurious contribu- 
tion from the CMB. 

The peculiar radial velocity of a cluster is given by 

m e c 2 y 

The relative error in the velocity is thus expressed as fol- 
lows: 

Sv r 5 A 5T e 

where A = AT / T . The error due to the CMB contam- 

v 

ination appears in the measurement of the temperature 
fluctuation AT/T and it is measured using a spatial fil- 
ter that optimizes the signal to "noise" ratio (the main 
"noise" being the primary CMB). The AT e /T e term is 
the error due to the uncertainty on the intracluster gas 
temperature determination. It should be derived from X- 
ray data. Hereafter, we do not include it in the evaluation 
of the ^v. 

To minimize the contamination of the CMB, we con- 
struct a spatial filter over which we calculate a variation 
of both AT/T and y on a single cluster. The spatial filter 
used hereafter computes the difference between the mean 
values of AT/T and y in the central part of the cluster 
and their mean values taken in a ring around the peak. 
This filter has three free parameters, one is the radius of 



the central disc and the two others the inner and outer 
ring radii. 

If the noise due to the spurious signals is dominated by 
a single component of a known power spectrum, the spa- 
tial filter can be optimized (Haenhelt & Tegmark 1996). 
In practice, several sources contribute to the noise includ- 
ing some non gaussian ones like the confusion of weaker 
clusters. Furthermore, we adress the problem of the opti- 
mization of the filter using only the informations included 
in the data. The shape and size of the optimum filter de- 
pend on the extent of the cluster and its density profile. 
We thus empirically optimize the simple filter in the case 
of one cluster model (King profile) and one primordial 
CMB spectrum (standard CDM with fl b = 0.05, h = 0.5 
and flo = 1) to evaluate the overall accuracy of the veloc- 
ity determination. The errors are dominated by the modes 
of the CMB power spectrum of wavelengths comparable 
or greater than the core radius. 

The optimum parameters of the filter are obtained by 
minimizing 5v rms , where the rms velocity is obtained af- 
ter many realizations of the millimeter and sub-millimeter 
sky including the astrophysical contributions due to CMB, 
foregrounds and clusters of galaxies and the instrumental 
noise. 

5.2. Results 

Our goal is to have a geometrical filter that could be ap- 
plied for a wide range of cluster sizes and derived directly 
from observations. We find that the optimized spatial filter 
has a central disc, associated with the peak of the cluster, 
corresponding to the region where y is greater than 70% 
of its maximum value j/o- For the ring, the best compro- 
mize is obtained for an inner radius Ri n = 0.5 x FWHM 
and a width AR = 2 pixels. Hereafter, we use the same 
optimized filter parameters for all cluster sizes. We check 
that varying both AR, in a range of 1 to 3 pixels, and i?, n , 
in range of 0.5 to 1.5 x FWHM, introduces an error of a 
few percent only in the velocity determination. 
Figure 6 displays the Sv rms obtained with the optimized 
filter (for a cluster with yo ~ 10~ 4 ) as a function of the 
core radius of the cluster. It shows that the rms velocity 
decreases for decreasing sizes of the clusters. We expect 
that the modes of the CMB power spectrum which con- 
tribute mostly to the measurement of the radial velocity of 
the clusters are those with wavelengths comparable to the 
size of the ring of the geometrical filter, which is given by 
2 x (R in + AR) — (3.5 x 9 C + 6'). Thus, the contamination 
of the CMB increases with the size, as more contribution 
from the first Doppler is included. This explains the rise 
of Sv rms from 2 — 3 arcminutes upwards. For clusters with 
core radii smaller than 2 arcminutes, the beam dilution 
leads to a fast degradation of the velocity determination 
accuracy. We use gaussian beams of 9 and 5 arcminutes 
in order to evaluate the effect of the beam dilution. As 
expected, Figure G shows that the accuracy improves for 
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Fig. 6. Sv rms due to the contamination by CMB fluctuations 
as a function of the core radius of the cluster, for a cluster with 
y w 1CP 4 . Solid line: convolution with a 9 arcminutes beam; 
dashed line: for a 5 arcminutes beam. 



a smaller beam, for all the cluster sizes but the effect is 
striking for core radii smaller than 1.5 arcminute. 

We check that the velocity determination accuracy for 
the largest clusters does not vary when we increase the 
noise since, in this range, the CMB contamination is dom- 
inant. For small clusters, we also find that the Sv rms docs 
not increase by more than 30% when we increase the noise 
by a factor 3. 

We compare the accuracy determination for a rich clus- 
ter (y — 10~ 4 ) and a weak one (y — 3. 10~ 5 ). Obvi- 
ously, the evaluation of the velocity is more accurate for 
the richest clusters. In fact, we find that, as expected, the 
velocity uncertainty scales as 1/yo- This enables to derive 
the accuracy in the velocity determination for each clus- 
ter, knowing its core radius 9 C and its central y parameter, 
t/o- In the framework of a particular cluster model (here 
Bartlett & Silk 1994), one can deduce 5v rms as a function 
of the cluster mass. 

Hacnhelt & Tcgmark (1996) have discussed in de- 
tail the optimized spatial filter for CMB contamination 
only, for various cosmological parameters. Our simula- 
tions, which take into account a more realistic noise to- 
gether with the limitations due to other astrophysical con- 
tributions (dust, free-free, ...), confirm that the accuracy 
they found can be achieved when the other sources of noise 
are taken into account. 

An additional uncertainty, noted 5T e /T e in Eq. 5, 
comes from the determination of the temperature of the 
intraclustcr medium from X-ray observations. This uncer- 
tainty is strongly dependent on the observed cluster and 
the characteristics of the instruments. With an instrument 
such as XMM-EPIC, the temperature of nearby clusters 
will be determined with a very high accuracy (~ 5%) for 
more distant clusters the uncertainty could be as high as 



10%. The results of the uncertainty on peculiar velocity 
evaluation, given Fig. 6, have been obtained neglecting 
the contribution of the uncertainty on the temperature 
T e . This contribution, which is not dominant, should be 
added according to the accuracy of the specific X-ray data 
used. 



6. Velocity dispersion measurements 

Our simulations show that the measurement of the pe- 
culiar cluster velocity is marginally possible only for the 
strongest clusters, in terms of their y parameter (y > 
10~ 4 ). Furthermore, the number of such clusters over the 
sky is too small to give very useful statistical information. 
Meaningful measurements can only come from a statisti- 
cal analysis, one beeing the rms velocity dispersion, the 
other beeing the bulk velocity. 

The measurement of the rms radial velocity 
(vrms) measured of a cluster is the combination of its real ve- 
locity dispersion v rms and of the velocity uncertainty due 
to the background contribution (v rms ) inst r- We therefore 
have: 



rms /measured 



"i" {Vrms) 



rmsJinsLr 



The error in the determination of the rms velocity on N 
clusters decreases like 1 / y/N. Therefore, it is necessary to 
measure the SZ effects on a large number of clusters in 
order to have the smallest error in {v rms ) mea8ure d- Never- 
theless, this kind of information is only partially relevant 
since the main difficulty is that such a measurement re- 
quires a good evaluation of (v rms ) instr which is not easily 
achieved. 

Over large scales, another accessible piece of statisti- 
cal information from the peculiar velocity of the clusters 
is the bulk velocity averaged over given volumes. In the 
specific model of this paper: De Luca et al. (1995) for the 
cluster counts and Bartlett & Silk (1994) for the evolution, 
and taking into account the limits due to the sensitivity 
of the experiment, one can derive the number of observed 
clusters per unit solid angle and intervals of redshift and 
mass. For a given mass and redshift the observable pa- 
rameters (9 C , y n ) are kown and using the results of section 
5.2. (Fig. 6), one can evaluate the accuracy 8v rms of the 
radial peculiar velocity for each class of clusters. 
In a given volume containing N clusters with individual 
peculiar velocities Vi and accuracies ct^, the best estimate 
of the bulk velocity V is given by the mean weighted ve- 
locity: 



N 1 



V =< v >- 



i=0 



N 1 



i=0 
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One can compute the overall accuracy a in the same vol- 
ume, which is given by: 

2 t=0 1 

We define a local volume V; oc by the volume being within 
the redshift range < z < 0.05 and which corresponds 
to a radius of 150/i _1 Mpc. We also define, up to redshift 
z = 0.9, volumes equal to V (oc (V; oc ~ 10 7 /i~ 3 Mpc 3 ). 
We compute the number of clusters in these volumes by 
integrating the counts over masses and redshifts. For each 
cluster, we also compute the observables C and y . Using 
the results given Fig. 6, to which we add a conservative 
uncertainty of 20% due to the estimate of the intracluster 
temperature, and together with the computed values of 6 C 
and yo, we derive the individual accuracy in the pecliar 
velocity determination for each cluster. We thus compute 
the overall accuracy a in each volume. 

Our estimates are summarized in Table 1. We find that 
the overall accuracy in the local volume is about 60 km/s. 
It goes through a minimum around z = 0.1 but it lower 
than 100 km/s for 0.5 < z < 0.7. At higher redshifts, 
the overall accuracy in the peculiar velocity determination 
reaches about 200 km/s. 



Table 1. Overall accuracy a in the peculiar velocity measure- 
ment for different redshift ranges. 



z 


fi (sr) 


a (km/s) 


0-0.05 


47T 


60 


0.05-0.1 


2.2 


12 


0.1-0.3 


0.157 


28 


0.3-0.5 


0.098 


34 


0.5-0.7 


0.096 


94 


0.7-0.9 


0.11 


223 



The measurements of the the bulk velocities from the 
peculiar velocities of individual clusters obtained using the 
combination of the SZ thermal and kinetic effects can be 
done with an overall accuracy better than 100 km/s up 
to z = 0.7. We have thus shown that a survey of the sky 
in the SZ effects followed by a redshift survey of a large 
number of detected clusters gives the possibility of map- 
ping the velocity fields in the Universe on very large scales 
(> lOO/i -1 Mpc) with a good accuracy. In the assumption 
of a potential flow and using some reconstruction method 
such as POTENT (Bcrtschinger & Dekel 1989), one can 
then derive the full three-dimensionnal velocity field and 
thus the density field which in turn allows comparisons 



with the one traced by galaxies. Therefore, the applica- 
tion of the SZ effect measurements to the evaluation of 
the peculiar velocity of clusters of galaxies should become 
a very useful tool to test and constrain theories of struc- 
ture formation and evolution. 

7. Conclusions 

Future high sensitivity and high resolution CMB ex- 
periments will have to remove the contributions from 
various foregrounds, including that generated through 
the Sunyaev-Zeldovich effect. Fortunately this effect can 
be easily separated through its spectral signature. Con- 
versely, this offers the exciting prospect of creating a 
large catalogue of clusters selected entirely via this effect, 
thereby allowing detections at high redshifts. In addition, 
using both thermal and kinetic SZ effects, one should then 
be able to map the very large scale velocity field. 

In order to make rather realistic predictions on the 
potential of such experiments, we have used complete 
sky simulations developped to assess the capabilities of a 
multi-frequency (30 to 800 GHz) space mission dedicated 
to the observation of the CMB between a few arcminutes 
to ten degrees, with a sensitivity close to AT/T ~ 2. 10 -6 
(Bouchct et al. 1995). These simulations take into ac- 
count both expected astrophysical components and the 
instrumental characteristics of the future space mission 
COBRAS/SAMBA. 

We have shown that a COBRAS/SAMBA like experi- 
ment will indeed give a fairly complete catalogue of more 
than 10 4 resolved clusters, up to a redshift of one or more. 
In the case of resolved clusters with central comptonisa- 
tion parameter y > 2. 10~ 6 , it is possible to reconstruct 
their y profiles up to one degree radius for the strongest 
ones, showing that the main limitation in the observation 
of the outer parts of the cluster profile will be due to the 
confusion with weaker clusters. 

Measuring the peculiar velocities of clusters is possi- 
ble (Eq. 3) when we combine both SZ thermal and kinetic 
effects. We have constructed and optimized a geometrical 
filter for this purpose. Our results, taking into account 
all the astrophysical and instrumental contaminants in 
addition to the CMB emission, confirm the Hcanhelt & 
Tegmark (1996) results for which CMB was the only spu- 
rious signal. We have also shown that measuring the pecu- 
liar velocities, for a large number of clusters up to z = 0.5 
gives overall accuracies on bulk velocities, computed in re- 
gions of typical dimension 100 hr 1 Mpc, better than 100 
km/s. This suggests that it will therefore be possible to 
map the density field of the Universe using SZ velocity 
measurement. 

The overall accuracy at higher z depends more strongly 
on the cluster counts which are not constrained strongly. 
In any case, the specific model we used here, which takes 
into account the evolution of the number of clusters gives, 
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if anything, an underestimate of the counts compared with 
no-evolution models. 



We have seen in, Sect. 2.1, the relation between the y 
parameter and temperature variation which gives: 
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8. APPENDIX A: y profile for a resolved cluster 

In the Rayleigh- Jeans part of the spectrum, the SZ ther- 
mal effect is given by: 



1 1 AT 



ATrj = -2T C 



MB- 



kUT 

m e c 2 



T e dl. 



Using the assumptions of isothermality and spherical sym- 
metry to describe the gas distribution, we have: 

ATrj = -2T C MB ^n e0 T ea DA f fnd(, 

m e c 2 J 

with: 



fn = 1 + 



2 i /*2\ - — 



+ c 



6 2 



We have R m ax = 20 R c which we take as an integration 
limit for the integral above. 

To compute the integral up to R m ax, we split the in- 
tegration to infinity into two parts and write: 



Io,R — Ia,oo — Ir,, 



with: 



q2 ?2 \ -« 



Ir,oo — 



1 



d a jr„. 
We finally find: 



i 92 l2 Y\n 



j - l o /- r(t*-i/2) / e 2 



()2\ , ,-(l+2a) 
u \ u-> max 



1/2- 



,l-2a 



1 -2a 



2 V 1 + 6*2 y l + 2a 



where: 



r>2 

2 -1 1 lx max 

max ' jj>2 



Therefore, the temperature variation ATrj is written as: 



kT, 



ATrj = —2TcMB ^(TTrieoD Ah.R- 

m e c z 



y 



RJ 



2 fix) T C MB ' 



The profile of the y parameter is thus derived directly 
from the expression of the temperature variation, and the 
AT/T profile is obtained using Eq. 2, Sect. 2.2. 
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